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Abstract 

The ti system allows a truly perturbative treatment of the potential. Completing 
previous computations, we calculate the contributions of QCD box graph correc- 
tions, which make the "relativistic" 0(a 4 ) corrections in the non-Abelian case differ 
from the well known corrections of the same order in QED. 
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1 Introduction 



The literature on the perturbative treatment of the quark-antiquark potential is very 
abundant. But even for the bottom quark these calculations suffer from the drawback that 
non-perturbative contributions can be estimated to be as large as the leading Coulomb 
part, and thus that their perturbative inclusion is not justified. On the other hand, 
phenomenological potentials with a large number of free parameters are used extensively 
0. The situation is completely different for the top quark. With its mass of ~ 180 
GeV it is heavy enough to allow a completely perturbative treatment. The drawback 
of this heavy mass is that, owing to its weak decay, individual bound state levels are 
not directly observable. However, this turns into a virtue since this large decay width 
effectively cuts off non-perturbative contributions in the calculation of the cross sections 
near threshold ||. For computations in this region the input of a perturbative QCD 
potential is justified and it is thus desirable to investigate to which extent we are able 
to calculate such a potential. Unfortunately, despite the large amount of literature, the 
number of studies including reliable QCD perturbative calculations beyond the lowest 
orders is actually small. For example in most cases higher order corrections are derived 
from on shell scattering amplitudes 0. This is not a consistent approach since different 
powers give contributions to the same order in the energy shift of bound states. The 
present article is also intended to shed some light on the calculation of spin-dependent 
forces up to 0(a^ In oc B ) where perturbation theory is likely to be applicable even for lighter 
quarks. 

Recently a consistent treatment on the basis of the Bethe-Salpeter (BS) equation 
was given where the level shifts themselves have been used as ordering parameter. 
There it was emphasized that it would be incorrect to extrapolate from the absence of 
corrections to 0(af) other than the tree graphs in the QED case to QCD because gluon 
splitting allows new types of graphs. Therefore we will discuss here possible sources of 
new corrections to that particular 0(af) in some detail. 

2 QCD two-loop box graphs 

The two-loop box graphs on which we concentrate in this letter are represented in fig. [I]. 
As we adopt the Coulomb gauge, the full gluon propagator (graphically represented by 
a wavy line) can be split, as usual, into a Coulomb part (broken line) and a transverse 
part (curly line). Contributions to the two-loop box graphs, which do not appear in fig. 
[3], may be checked to give higher-order corrections. 

In H|] an exact zeroth-order solution has been derived for the toponium BS equation, 
generalizing the Barbieri-Remiddi equation || to the case of unstable particles. For our 
purposes it is sufficient to consider the corresponding zeroth-order wave functions x an d 
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X in a non-relativistic approximation: 



X (p)» = _^!L^ m (p) 5 = - 7 oX(p) nr 7o , (2.1) 

where uo n := (p 2 /m — E n )/2 — ie, E n = —ma 2 /4n 2 are the Bohr levels, a := (4/3)g 2 /47r = 
(4/3)a s , n ; m the Schrodinger-Coulomb wave functions, x an d X are written as spinor 
matrices with S a 4 x 4 matrix representing the spin state of the particle-antiparticle 
system: 




singlet 

triplet (m = 0, ±1 a? m = 1) 



S = 

X ± := (2.2) 



All the following contributions are evaluated between S and — 70S70 (times a normaliza- 
tion factor 1/6). 

First we will consider the QCD box graph of fig. |l|a. Let pi, —P2, be the incoming 
and p[, —p' 2 , the outgoing momenta of the particle and antiparticle. We define the total 
momentum P := pi — p 2 = p[ — p' 2 , the relative momentum of the incoming particles 
p := (pi + P2V2 and q := p' 2 — p 2 = p[ — p\. Our present correction is put between BS 
wave functions, in agreement with the general approach of BS perturbation theory. In 
that formalism, after performing the colour and spinor indices traces, the perturbation 
kernel of graph |T]a can be written as 



= 12 ig & J 



d t d 4 k (p\ — to + m) (p 2 — k — m) 
(2tt) 4 (2vr) 4 [( Pl - t) 2 - m?][(p 2 - k) 2 - m 2 ] 

' Q(q-k,k,t-k), (2.3) 



^t-k) 2 k 2 (q-k) 2 t 2 (q-t) 



with 



Q(p, q ,k): = pq- ^ ^ tL . (2.4) 



(p-k)(q-k) 
k 2 

In the centre-of-mass frame P = 0, p\ = (Pq/2 +Po,p) and p 2 = (—P /2 +Po,p)- The 
integrations over t and k are quite simple, but the complicate angular structure seems 
to be a serious obstacle to the exact evaluation of the remaining integral. Observing that 
for the bound state the scales of the momenta [] are p, q ~ ma and po ~ ma 2 , it is natural 

1 All calculations have been performed with the help of the program of symbolic manipulations FORM 

i- 

2 This is the reason why a simple counting of vertices fails when estimating leading-order contributions 
of the Feynman graphs to the levels of bound states. 
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to introduce rescalings as follows: 



Po 
P 

q 
t 
k 



a 2 p , 
ap, 
aq, 
a t , 
a k . 



(2.5) 



As a consequence, the right-hand side of eq. (|2.3j ) is easily seen to retain a single power of 
a times an integral which is finite in the limit a — > 0, remembering that the centre-of-mass 
energy is Po — 2m + 0(a 2 ). Therefore, the leading contribution in a of ( |2.3|) is simply 
obtained by putting a = in this integral. The resulting expression is 



162 7r a a 



dH d 3 k 
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(2vr) 3 (2tt) 3 fc 2 (q - k) 2 t 2 (q - t) 2 (t - k) 2 



Q(q-k,k,t-k) (2.6) 



which only depends on q. The trick to simplify this expression is to write the inner 
products in terms of quadratic expressions that partly cancel the denominators: 



Q(q-k,k,t-k) 



2 q 2 - k 2 - t 2 - (q - k) 2 + (k - t) 2 - (q - t) 



t 2 (q-k) 2 t 2 (q-t) 2 A; 2 (q-k) 2 _ fc 2 (q-t 
(k-t) 2 + (k-t) 2 + (k-t) 2 (k - 1) : 



(2.7) 



While the original integral (|2.6| ) is finite, the integrals one obtains from a single summand 
in ( |2.7|) are not, so the evaluation of them requires some care. In the appendix we 
give the list of integrals necessary to evaluate (|2.6|) . The divergent integrals have been 
regularized with a cut-off in the integration region, but, as a check, we have verified 
that the same result can be obtained within dimensional regularization. Removing the 
auxiliary regulators, obviously all divergent terms cancel in the final expression, which 
can be written as 



81 , 

7T 12 

128 v 



7T 



a" 



(2.8) 



where we have restored the initial definition of q (before scaling). Taking into account 
the factor a 3 from the wave functions, we find, as a consequence of eq. ( |2.8|) , that the 
graph [T|a indeed yields contributions of 0(a A ) to the energy levels. Therefore it must be 
included in an 0(a 4 ) evaluation of the QCD potential V. Actually our result confirms 
some previous qualitative estimates 0-[f|. 

Graph |I]b produces contributions to the energy levels that are of higher order than 
|l]a. In fact the insertion of a transverse gluon propagator between the wave function and 
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Figure 1: Two-loop box graphs. 



the free quark propagator S annihilates its leading contribution. But additional powers 
in the spatial momenta and consequently (see ( |2.5|) ) in powers of a remain: 

s (4 + ") - ± p 0+ P /L Ep+ 4 A± -^ + -) + -' < 2 - 9 > 

( ± f +P ) « ^ i ±po+P J 2 _ Er+i£ (-^) + ■ ■ ■ ■ (2-10) 

with := i/p 2 + m 2 . However, we note that this graph contributes to the spin- 
dependent part of the interaction. 

Graph [TJc was considered in |JJ] , where it has been shown that the corresponding leading 
contribution to the energy level shift is 0(a 5 \na). 

Graphs like fig. |l|a and |l|b with crossed Coulomb lines (fig. [l|d) are irrelevant because 
they vanish due to group theoretical factors. 

By connecting the quark lines between the interactions in the non-Abelian vertex 
correction graph and graph |l|a, as shown in fig. [| one obtains graphs which could be 
named non-Abelian Lamb shift graphs. Consider for instance the graph of fig. 0a with 
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only one Coulomb gluon insertion, and call t the momentum carried by the transverse 
gluon, s and k the momenta of the Coulomb gluon connecting it with the upper fermion 
line. Taking into account only the leading contribution to the fermion propagator (the 
first term of eq. Q2.9Q ), we obtain 

dH d 3 k d 3 s 1 



= -81ttV 



(2tt) 3 (2tt) 3 (2tt) 3 t(s - t) 2 F(k + t) 2 s 2 (q - k - s) 2 



(t + E p ik + E p >k t — P ) (t + E ps + E pst — P ) ® ' ' )' ( ) 



with E p , k := yj (p' - k) 2 + m 2 , E v < u := ^(p 1 - k - t) 2 + m 2 , E ps := yj (p + s) 2 + m 2 and 
Ep S t '■= \J (p + s — t) 2 + m 2 . Scaling all momenta according to eq. (|2.5| ), Hell takes the 



form 

K 1 ) _ qi^4_2 



^2.a = -8l7rVJ(a,q), (2.12) 



where / is the integral of eq. ( 2.11 ) after scaling. Since it can be shown that the leading 
contribution of the derivative of I with respect to a is proportional to 1/a (scaling t a 
second time and putting a = in the integral), the leading contribution of H^l is 

= tg a o 2 lna. (2.13) 

From this argument it follows that all the graphs of fig. will contribute at least at the 
usual Lamb shift level 0(a 5 lno;). 




Figure 2: Non-Abelian Lamb shift graphs. 



3 QCD potential and renormalization schemes 

The 0(a 4 ) correction from the non-Abelian box graph fig. [TJa leads to a slight enhance- 
ment of the attractive Coulomb force. This is the only contribution of this type in our 
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renormalization scheme, which is defined as follows: m t is the pole mass, Z g i uon and 
Z\f, coulomb are defined by a subtraction at \i and 0, respectively. The updated potential 
given in M now reads 



V 



4m 3 m 2 



" / p2 ! r l r PJP 



2m 2 r 



3a 



2m 2 r 3 
33a 2 



LS + 



2m 2 r 3 



87rr 
81 



(7 + ln//r) + 



a 



128 



vr 12 



a' 



' 3(rS) s 
2 5 

E 

2a 3 



+ 



47ra 
3m 2 



S 2 5(r) 



47rr 



Ei( 



-rrrijes 



5 1, 
- + - In 

6 2 



u 5 



+ 



9a 2 



7T~ 



47rr (167r) 2 r 



27 2 



7T 



^- + 2(7 + In fir) 2 

D 



576(7 + l n 



8 47raQED(/i)a 
"9 r 



^2 GFm 2 ^! + v^G F m|a 2 ^ f 7 - -S 2 ) 
47rr J m| V 3 / 



+v / 2G F m|( 



2nr 



:sr) 



2a 2 



rrizr 



+ 



m 7 r 



:sr) 



(3.1) 



where m^, m# and are the quark, Higgs and Z masses. 

However, the usual predictions for scattering measurements are in the MS scheme. 
Thus a short remark on the coupling constant to be used seems in order. It is customary 
to take the strong coupling constant a s in the MS scheme. In particular experimental 
determinations are always given in terms of a s MS (mz) [10]. But in bound-state calculations 
it is natural to use an a s defined differently. This has been done in Qj. Clearly it should 
be possible to relate the two schemes. 

The quark-antiquark potential in the MS scheme with only light quarks to 0(a s ) in 
the self-energy function of the Coulomb gluon can be taken for instance from [p][ . The 
contribution of heavy quarks in this scheme has been calculated in [12|. We can thus 
write for the potential in momentum space 



47ra 



33, 31 



47T \ 4 



1=1 



2 ' 

m z A 5 



6 2 \fi 2 n 



(3.2) 



Comparison with the potential V gives 



a 



BS 



a 



MS 



Y _ ft MS 
167T 



31-E 



i=l 



10 



- 2 In 1 + 



/// 



// 



'3 - 



(3.3) 



where BS denotes our bound-state scheme. It should be emphasized that we have explic- 
itly included the dependence on massive flavours in this formula. Numerically, eq. (|3.3|) 
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means that our a is slightly smaller than the usual one, which is advantageous for our 
perturbative calculation. For m t = 180 GeV, a^ s (m z ) = 0.119 ± 0.006, we obtain for a 
renormalization at /i = (re = 2 /ma is the Bohr radius of the bound state): 

/I = 17.5 ±1 GeV, (3.4) 
a m (p) = 0.205 ±0.01, (3.5) 
a(fi) = 0.19 ±0.01. (3.6) 

Remember that the above values for a, a MS differ by a factor of 4/3 from a B . 

For the comparison of the two schemes, the terms of the form a™/q 2 have been im- 
portant. To compare to relative 0(a 2 ), terms up to af/q 2 are needed in the potential. 
While the only term of this form in our scheme has been calculated in our present work 
(eq.( ^.8p ), we are not aware of an analogous calculation of the terms of order a 2 /q 2 in 
the MS scheme. Therefore, at present the strong coupling constant in these two schemes 
can only be related to 0(a s ). 

4 Conclusions 

In this note we present an analytical calculation of the contribution of the QCD box 
graphs to the potential that would give rise to an 0(a 4 ) energy level shift (eq. ( |2.8|) ). 
This result is new and confirms previous qualitative estimates (0-||). It provides an 
example of a non-trivial contribution to 0(a 4 ) of a two-loop graph, which certainly is 
not included e.g. in the running coupling constant and which is typical for a non-Abelian 
theory. 

While the calculated correction appears in our bound-state scheme as well as in the 
more common MS scheme, some corrections of the same form have not been calculated 
as yet to the required order in the latter. Thus a direct comparison of the two schemes is 
out of reach, at present. 

We also identified many non-Abelian graphs that give rise to 0(a 5 In a) contributions. 
These involve two-loop graphs as well as infinite chains of graphs, which could be called 
"non-Abelian Lamb shift graphs" . 
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gramme "Human Capital and Mobility", Network "Physics at High Energy Colliders", 
contract CHRX-CT93-0357 (DG 12 COMA). 
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A Appendix 



In this appendix we give some of the integrals needed to obtain ( |2.8|) from ( [2.6|) . The di- 
vergent integrals are regularized with cut-offs in the integration region. For the evaluation 
of these integrals within dimensional regularization we refer the reader to Jl3[ : 

dH 1 1 1 111 



(27r) 3 t 2 (q-t) 2 (t-k) 2 8g£;(q-k) 2 ' 
dH 1 1 11 



(A.l) 
(A.2) 



(2tt) 3 (q- t) 2 (t -k) 2 8|q-k|' 
r d 3 k 1 1 11 

L = Sv (2hlfe) " " l(fl) " lnfe) " ' (A - 3) 

L~ wr^i^w = i> (?) - ln(£2)) • (A - 4) 

L wf^w^w = i? (ln(,) " ln(ei) + 11 • (A ' 5) 

/ (fliPi"^ = AV(9) - ln(ei) + 1) , (A.7) 
Jk>t 1 (2tt) k |q — k| 271^ q 

( d 3 k 1 1 11 

L, (2^*3 |^kp = W< 2 "^ " ^ " h < £2 » " ^ ' < A - 8 > 
Finally a useful identity is 

d 3 A; dH 1 / 1 1 \ 1 11 



(2tt) 3 (2tt) 3 k 2 \{q- k) 2 (q - t) 2 J (t - k) 4 32vr 2 g 2 ' 



(A.9) 
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